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Abstract. We present a method of mode analysis to search for signals with frequency 
evolution and limited duration in a given data stream. Our method is a natural 
expansion of Fourier analysis, and we can obtain the information about frequency 
evolution with high frequency precision and high time resolution. Applications of 
this method to the analysis of in-spiral and burst signals show that the signals are 
characterized by an index which we name "bounce hardness" . The index corresponds 
to the growth rate of the signals. 
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1. Introduction 

Interferometer gravitational wave detectors have been operated around the world to open 
a new window for astronomy [TJ [2j El IH El [7] . Though gravitational waves have still not 
been detected, upper limits for the wave signals such that from binary coalescence [H HI |5] 
and from stellar core collapse [HI |7] have been obtained by the operations. 

To search for the signals, matched filtering, which is a method with predicted 
waveforms, is commonly used. For example, exponential Gaussian correlator [HI |9] uses 
Gaussian envelope, O-pipelinepO] uses bi-square envelope, and coherent WaveBurst [TTj 
uses Mayer wavelets. Some of them are used with a variation of linear predictive 
coding[12], which was originally developed for speech recognition. 

Matched filtering is successfully used to analyze not only the chirp signals by binary 
coalescence but also the burst signals by core collapse [8J, where the waveforms caused 
by in-spiral binary stars are well anticipated by post-Newtonian calculations, and where 
there are various waveforms caused by stellar core collapse [13] depending on the models 
of gravitational wave emission process. 

In near future, the signals will be found by the matched filtering approach. Once 
the signals are obtained with certain signal-to-noise ratio (SNR), a method for detailed 
analysis of the signals is required. Obtaining the characteristics of the observed waves 
by the detailed analysis, more efficient signal detection comes to possible. 

Therefore, we introduce our method of signal analysis for the detailed analysis. 
Our method is a kind of mode analysis, in which we fit a given time series with a linear 
combination of nonlinear oscillators. Our method with certain limitations is equivalent 
to Fourier analysis, and without the limitations, the frequency precision and the time 
resolution of the frequencies are beyond those of Fourier analysis. For example, we can 
obtain the frequency information even from a fraction of cycle data, and we can obtain 
the characteristics of the data without oscillation, as we show in the following sections. 
As our method does not require specific assumptions such like predicted waveforms, it 
is applicable to the general field of signal analysis. 

Through the analysis, we show that the waveforms of pre-bounce gravitational 
wave signals are characterized and categorized by an index J, which we name "bounce 
hardness" . The index corresponds to the growth rate of the signals. 

To begin with, we introduce our method of analysis and its simple implementation. 
Then, we apply the method to the analysis of in-spiral and burst signals without 
background noise, and we show that the signals are characterized by an index of bounce 
hardness. Finally, we summarize the paper. 

2. Method of analysis 

The method we introduce here is a kind of mode analysis. 

Suppose that g(t) is the output of a dynamical system which we want to analyze, 
and suppose that we can obtain g(t) only as a numerical time series data. Then, we fit 
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the time series g(t) with a linear combination of oscillators with nonlinear oscillation. 
The fitting function h{t) is written as 

M 

h(t) = eX P H m(t) (1) 

m=l 

where M is the number of oscillators and H m (t) is the complex function of the nonlinear 
oscillation which is continuous and differentiable. 

For a short enough time scale, the fitting function h(t) is approximated by a linear 
function. The linear function ho(t) is written as 

M 

ho{t) = E ex P [H m {to) + H' m (to) (t - t )] (2) 

m=l 

for t ~ t - H m (to) and H' m {to) correspond to the power/phase and the mode of the 
oscillation, respectively. 

As the approximated function ho(t) is a linear function, we can easily fit the given 
time series g(t) with the function for t ~ to- 

We show an example of the fitting: in general, the time series g(t) is given as a 
discrete time series g(to + nAT), where n is the index of the time series and AT is 
the time interval of the sampling; we fit the discrete time series with a linear difference 
equation 

M 

exp [H' m (t )AT] D} g(t ) = (3) 

m=l 

where D is the time shift operator such that 

D n g(t ) = g(t - nAT) (4) 

and the mode of the oscillation H' m (to) is obtained by solving the linear difference 
equation (jHJ). 

We introduce the following functions 

c m (t ) = H m (t ) (5) 

Uto) = ^H' m (t ) (6) 

X m (to) = KH'Jto) (7) 

to give physical meaning that the power/phase c m (to), the frequency f m (t Q ), and the 
decay/growth rate \ m (t ) at time t = to. 

By use of the functions, the fitted function h(t) for t ~ to is written as 

M*)Lt„ - Mi) (8) 

M 

= J2 ex P { c m(to) + [i27r/ m (to) + A m (t )] (t - t )} (9) 

m=l 

and the time series of the oscillating mode is written as 

H'Jt) =i2n f m (t) + X m (t). (10) 
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We note here the relation with discrete Fourier transform (DFT). The following 
settings 

C m {t)=C m (11) 
777 

fm(t) = 7 — (12) 

A m (i) = (13) 
M = N (14) 

where N is the number of samples for the analysis, gives DFT 



h (t) = 2 ex P 



N n 

C n + i2n— — (t - t ) 
NAT 



n=l 



(15) 



itself. Therefore, it is found that our method of analysis includes traditional Fourier 
analysis. 

The major characteristic of our method of analysis is that we can obtain the time 
series of the oscillating modes ( FlOl) of given time series with high time resolution. 
Each mode contains the information about frequency and decay/growth rate. In 
contrast, Fourier analysis is a method which has fixed modes and fitting parameter C n . 
Therefore, time/frequency resolution is limited by the number of samples for analysis, 
and decay/growth rates are not obtained directly. 

3. Implementation 

In this section, we introduce a simple implementation of our method of analysis to 
apply the method to the analysis of gravitational wave signals. The precision and time 
resolution of our method of analysis are shown through the application. 

We employ least-mean-square (LMS) fitting for the analysis in this paper. The 
simultaneous equations for the fitting are 







o 

uu m > n=1 _ N 



M 



1 - ]T a m D m g(t + nAT) 



2 



(16) 



m=l / 

for m' = 1,2, ...,M, where a m is the coefficient of the linear difference equation to be 
fitted, and iV is the number of samples for the LMS fitting. We obtain a m by solving 
the simultaneous equations. Then, the oscillating mode H' m (to) is obtained by the 
factorization 

M M 

1 - £ a rn D m = J] {1 " exp [H' m (t )AT] D} . (17) 

m=l m=l 

We solve the simultaneous equations for 

t n = t + nAT, (18) 

and we obtain the time series of the oscillating mode H' m (t n ). 

We note here that this implementation corresponds to the traditional maximum 
entropy method (MEM) when M and iV are large enough. Therefore, background 
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random noise, whose time series can not be written by a linear equation, is weeded out 
as the residual of the LMS fitting. In addition, we notice here that we should solve ffT6l) 
exactly. It is for the reason that the traditional method of the LMS fitting contains 
implicit periodic boundary condition to reduce computational complexity, and it gives 
wrong solution. 



4. Application to chirp signal 

In the following, we apply the above implementation to the analysis of chirp signal. 

We employ the chirp signal of in-spiral binary stars with mass 1AM & and the 
distance 10 kpc for the analysis [H]. The equations for the chirp frequency f c (t) and 
the chirp signal h c (t) are 

1 /c 3 \ 5/8 / 5 1 \ 3/8 

= * (256/iM^tbounce-tj (19) 



£5/3 



h c (t) = (vrM to J c (t)) 2/3 sin [27r/ c (t) (t - t honncc )} (20) 



c 4 r 



where M tot = 2.8M , /i = 0.7M Q , and r = 10 kpc. The time series for the analysis g(t) 
becomes 

g(t n ) = h c (nAT + t honQce ) (21) 

where AT is the time interval of the sampling, and we choose 0.05 ms as AT. We show 
the time series (12T|) in figure [TJ 




Figure 1. Time series of chirp signal (p?Tj) for the analysis. 



The bounce hardness index of gravitational waves 



6 



As the chirp signal h c (t) contains a single function for frequency, we choose 2 as 
the number of nonlinear oscillators M. Then, the chirp signal h c (t) is written as 

h c (t) = exp H+(t) + exp H-(t) (22) 

where the complex functions H±(t) of the nonlinear oscillators are written as 

H ± (t) = -- In (t bounce -t)± i2vr • 133.5 {t houncc - tf s — 43 =F |i. (23) 

In these conditions, the frequency time series f(t) of the nonlinear oscillators 
becomes 

f(t) = ^(H' + (t)-H'_(t)) (24) 

= 133.5 • - (^bounce — t) (25) 

and the growth rate time series A(i) of the nonlinear oscillators becomes 

(26) 

7 

(^bounce ^) • (^^O 

lb 

Here, the second term in (12"B1) is a correction term. The term is required because the 
complex functions H±(t) have strong nonlinearity, and locally linearized functions are 
different from the complex functions themselves. 

We choose 8 as the number for the LMS fitting N, and we plot the results of the 
analysis in figure [2] (chirp frequency f(t)) and in figure [3] (growth rate A(£)). As we 
employed the parameters M = 2 and N = 8, the number of samples for calculating f(t) 
and X(t) becomes 10 (M + N). In addition, as we employed 0.05 ms as AT, the time 
width of a single analysis becomes 0.5 ms. 

Figure [2] shows the time series of frequency with a semi-log scale. The solid line 
shows the theoretical time series of frequency ( |25l) . and the plots by circles show the 
calculated one from time series (l2~Tj) . Circles are plotted every 12 ms to maintain 
readability. The calculated time series of frequency agrees strongly with the theoretical 
one. The difference between them is less than 2% for t — bounce < —10, and is less than 
10% for t — bounce < -1- As we assume a linear equation for a duration of 0.5 ms, the 
difference becomes larger approaching to the bounce. 

Figure [3] shows the time series of growth rate with a log-log scale. The solid line 
shows the theoretical time series of growth rate (1271) . and the plots by plus signs show 
the calculated one from time series f l2T|) . The trend and the values of the calculated 
time series of the growth rate agree fairly well with the theoretical one. The oscillations 
observed in the time series of calculated values come from the LMS fitting. As our time 
width of analysis (0.5 ms) is much shorter than the cycle of oscillation of the chirp signal 
(~5 ms) in figure (TJ the fitted value is affected by the status of the curvature of the 
chirp signal. We can attenuate the amplitude of the oscillations by employing large 
enough iV compare to the period of the oscillations of the time series. However, as we 
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Figure 2. Semi-log plot of chirp frequency. Theoretically expected frequency (|25[> is 
shown by solid line, and calculated one from time series (j!?Tj) is plotted by circles. 
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Figure 3. Log-log plot of chirp growth rate. Theoretically expected growth rate (p?T|) 
is shown by solid line, and calculated one from time series (pHI) is plotted by plus signs. 
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assume a linear equation for a duration of (M + N) AT, it results low time resolution 
for tracing non-linearity. 

As shown above, our method of analysis gives information about frequency and 
growth rate with high precision and with high time resolution. 

5. Application to burst signals 

In this section, we employ the signals of gravitational waves calculated and presented 
by Dimmelmeier et al. [I5j [16] for the analysis. Among the signals, we employ the 
ones named A1B3G3 and A3B5G4. They represents two kinds of signal growth, and 
the other signals are categorized in either group. 

We prepared their time series data [T7j with 8 kHz sampling, and we employ the 
same implementation and the same parameters (M = 2 and N = 8) used in the previous 
section, because it is the most simple implementation. Then, we set the time series for 
the analysis g(t) as 

g(t) = A™ (28) 

which is given by Dimmelmeier et al., and we fit the time series with 

h{t) = exp H x (t) + exp H 2 (t). (29) 

As the pre-bounce time series has no oscillation, the imaginary part of H m (t) 
becomes zero, and we obtain two growth rates X\(t) and \2{t). 

Through the analysis, we found that their time series are characterized by an index 
of bounce hardness /, and the time series of growth rate H' m {t) is written as 

H'Jt) = X m(t) OC (bounce ~ • (30) 

Therefore, the fitted time series h(t) of gravitational wave becomes 

2 

h(t)= J2 ex P 

m=l 

where c m and d m are the fitting parameters for each time series, except when 1=1. 
The case 1 = 1 corresponds to the case in the previous section ( 1271) . where the real part 
of H m (t) is logarithmic. 

The signals with / < has no singularity at the bounce. In contrast, I > has 
singularity in the gradient, and / > 1 has singularity in the amplitude. In this way, the 
index / represents the hardness of the bounce. 

In the following, we show the results of the analyses for A1B3G3 and A3B5G4, the 
harder case and softer case, respectively. In addition, we categorize the other signals 
in the catalog [UJ by using the index /. 

5.1. A1B3G3 

Figure H] shows the time series of A1B3G3 and the fitted time series. Figure [5] shows 
the plots of growth rate A m (i) in a log-log scale. 



Cm ~\~ d m (tbounce £) 



(31) 
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Figure 4. Time series of A1B3G3 is shown by solid line, and our fitted time series 
(f3"4"| is plotted by circles. 
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Figure 5. Log-log plot of growth rate A. Two lines shown by plus signs imply 
polynomial functions. Two fitted functions (|3"1?]) and (|3Ul) are shown by dashed lines. 
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The plots of growth rates in figure [5] show two lines, and we fit them with straight 
lines. By the fitting, we obtain the following two functions of growth rates 

Ai(i) = 2.3(tbouncc _ t) °' 4 (32) 
A 2 (t) =0.25(t bounce -t)- 04 (33) 

and the fitted function in figure H] 

2 3 25 

h(t) = 702exp[- — (t bounce - t) - 6 ] + 168exp[-— r (t bouncc - t) a6 ] (34) 
0.6 0.6 

where we employed the obtained values at t = tbouncc — 0.4375 to set the parameters c m 
and d m in (I3T1) . The difference between A1B3G3 and (f34"|) is less than 10% for —1.4 < 
t — ^bounce < —0.1, and the maximum value of the difference for —10 < t — bounce < —1-4 

is 20% at t — ^bounce = —3.5. 

In this case, the index of bounce hardness becomes I = 0.4, and the gradient of the 
time series has singularity at the bounce. 



5.2. A3B5G4 

Figure [6] shows the time series of A3B5G4 and the fitted time series. Figure [7J shows 
the plots of growth rate A m (t) in a log- log scale. 

In contrast to the A1B3G3 case, the two lines plotted by plus signs in figure [7J are 
independent of time. This implies that X m (t) is a constant function. We fit them with 
straight lines, and we obtain growth rates 

Ai = 0.824 (35) 

A 2 = 0.108 (36) 

and the fitted function in figure [6] 

h{t) = 620 exp[0.108(i - bounce)] - 140 exp[0.824(t - £ bounce )] (37) 

where we employed the obtained values at t — bounce — 0.0625 to set the parameters 
c m and d m in (|3~T1) . The difference between A3B5G4 and fl37j) is less than 2% for 
—10 < t — bounce < 0, and is less than 4% for —22.5 < t — tbouncc < 1- 

In this case, the index of bounce hardness becomes 1 = 0, and the gradient of the 
time series has no singular points during the bounce. 



5.3. Categorizing signals 

We show the index / of each signal in the catalog[17J, and that of in-spiral binary stars 
in Table [U 

The table shows that the rapidly (large j3 TOt ; n i) and highly differentially (small A) 
rotating initial models exhibit exponential signal growth. Among them, the growth rate 
of A2B4G1 shows some fluctuation around a value, and it is not a simple constant. 
Therefore, it seems that a small enough A is required for stable exponential signal 
growth. 

Most of the other burst signals exhibit intermediate (between exponential and 
polynomial) signal growth, and the few resting signals are not identified simply. 
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Figure 6. Time series of A3B5G4 is shown by solid line, and our fitted time series 
(|5T|) is plotted by circles. 
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Figure 7. Log-log plot of growth rate A. Two lines shown by plus signs imply 
constant functions. Two fitted functions and are shown by dashed lines. 
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Table 1. Bounce hardness index / of signals. 





Signal growth 


Examples of signals 


/ = 


Exponential 


A2B4G1, A3B4G2, A3B5G4, A4B4G4, A4B4G5, A4B5G5 


< / < 1 


Intermediate 


A1B3G3, A1B3G5, and most of the other signals 


/ = 1 


Polynomial 


In-spiral binary stars 



6. Conclusions 

We introduced our own method of signal analysis, which is a natural expansion of 
Fourier analysis. Using our method, high resolution information about frequency and 
decay/growth rate is obtained from small number of samples. Therefore, the detailed 
analysis of frequency evolution and decay/growth rate evolution comes to possible. 

We applied the method to the analysis of the chirp signal and the burst signals of 
gravitational waves. We found that the signals are characterized by an index of bounce 
hardness /. The growth rate A m (t) of the signals is written as ( 1301) and the index 
/ is 1 for the chirp signal (polynomial signal growth), 0.4 for A1B3G3 (intermediate 
signal growth), and for A3B5G4 (exponential signal growth). The rapidly and highly 
differentially rotating initial models belong to A3B5G4 group, and most of the other 
burst signals belong to A1B3G3 group. 

We note here that this is the first paper introducing our method of analysis, and 
there are no references on the method. In addition, the method of analysis is still under 
development, and more qualified analysis will become possible through the progress of 
the development. For example, the way to detect signals from low SNR data stream by 
using our method of analysis, and the way to apply our method of analysis into matched 
filtering will be found. 

References 

[1] Kawamura S 2010 Class. Quantum Grav. 27 084001 

[2] Accadia T and Swinkels B L 2010 Class. Quantum Grav. 27 084002 

[3] Abbott B P et al 2009 Phys. Rev. D 79 022001 

[4] Abbott B P et al 2009 Phys. Rev. D 79 122001 

[5] Abbott B P et al 2009 Phys. Rev. D 80 047101 

[6] Abbott B P et al 2009 Phys. Rev. D 80 102001 

[7] Abadie J et al 2010 Phys. Rev. D 81 102001 

[8] Clapson A C et al 2008 Class. Quantum Grav. 25 035002 

[9] Acernese F et al 2009 Class. Quantum Grav. 26 085009 

[10] Chatterji S et al 2004 Class. Quantum Grav. 21 S1809 

[11] Klimenko S et al 2008 Class. Quantum Grav. 25 114029 

[12] Itakura F and Saito S 1970 J. IEICE 53A 35 

[13] Ott C D 2009 Class. Quantum Grav. 26 063001 

[14] Schutz B F 1986 Nature 323 310 

[15] Dimmelmeier H, Font J A and Miillcr E 2002 Astron. Astrophys. 388 917 



The bounce hardness index of gravitational waves 13 
[16] Dimmelmeier H, Font J A and Miiller E 2002 Astron. Astrophys. 393 523 

[17] URL http://www.mpa-garching.mpg.de/rel_hydro/wave_catalog.shtml MPA Garching Gravita- 
tional Wave Catalog 



